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ABSTRACT 

A  study  of  the  fluctuating  velocity  field  in  the 
developing  region  of  a  pulsating  pipe  flow  is  presented. 

Theoretically,  three  separate  determinations  of  radial 
and  axial  perturbation  velocity  profiles  were  calculated 
involving  varying  degrees  of  analytical  simplification. 

Experimentally,  pulsating  disturbances  were  created 
by  a  sinusoidal  pressure  wave  exciter  installed  at  the  out- 
let of  a  1.5  inch  circular  wind  tunnel  facility.   Cross- 
wire  anemometer  measurements  of  the  pulsating  velocity 
components  compared  favorably  with  the  parallel  theoretical 
results.   Likewise,  experimental  determinations  of  the 
phase  angles  between  velocity  and  pressure  signals  compared 
favorably  with  analytical  values. 
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NOMENCLATURE 


Constants,  coefficients  of  the  polynomial 
approximation  of  the  velocity  profile 


d  Inside  diameter  of  the  pipe 

P  Pressure 

r  Radial  coordinate 

rn  Radius  of  the  pipe 

R  Reynolds  number 

t  Time 

0  Mean  velocity  in  the  pipe 

UM  Maximum  velocity  in  the  pipe 

u  Axial  component  of  velocity 

v  Radial  component  of  velocity 

x  Axial  coordinate 

v  Kinematic  viscosity 

p  Density 

<J>  Angle  of  phase  shift 


ACKNOWLEDGEMENT 

The  author  wishes  to  express  his  profound  gratitude  to 
Professor  T.  Houlihan  for  his  guidance,  encouragement, 
and  infectious  enthusiasm  during  the  course  of  this  inves- 
tigation.  A  special  note  of  appreciation  is  also  given  to 
Messrs.  K.  Mothersell,  J.  McKay,   J.  Beck,  and  G.  Baxter 
of  the  Mechanical  Engineering  Machine  Shop  for  their 
generous  assistance  in  the  construction  of  the  experimental 
apparatus . 

The  author  is  especially  grateful  to  his  wife,  who  was 
patient  and  understanding  throughout  the  period  of  this 
study. 


I.   INTRODUCTION 

The  development  of  flow  in  the  inlet  region  of  a  circular 
cylindrical  tube  has  been  a  subject  of  interest  since  the 
time  of  Hagen  (1839)  and  Poiseuille  (1841) .   Osborne  Reynolds 
(1883)  conducted  extensive  investigations  into  the  stability 
of  pipe  flow,  and  this  area  has  become  a  classic  problem 
with  a  sizable  amount  of  literature  published.   However,  the 
study  of  fluid  flow  through  a  rigid  tube  under  the  influence 
of  a  periodic  pressure  difference  is  more  limited. 

The  theoretical  study  of  pulsating,  laminar  pipe  flow 
was  primarily  given  by  Th.  Sexl  [1]  and  S.  Uchida  [2]. 
Experimen^il  data  are  limited  to  the  works  of   Atabek,  Chang, 
and  Fingerson  [3]  and  the  results  of  Florio  and  Meuller  [4]. 
With  the  introduction  of  the   modern  computer,  several  attempts 
have  been  made  to  solve  the  problem  by  numerical  methods . 
The  most  recent  work  in  this  regard  is  that  of  Lew  and  Fung 
[5] .   In  the  theoretical  analysis  the  main  problem  has  been 
the  non-linear  convective  acceleration  term  in  the  Navier- 
Stokes  equation.   Uchida  [2]  was  able  to  apply  linear  theory 
for  the  fully  developed  laminar  pulsatinq  flow.   Lew  and 
Fung  [5]  found  this  linearization  to  be  valid  in  the  entry 
region  for  high  Reynolds  numbers.   However,  this  assumption 
leads  to  significant  errors  for  small  Reynolds  numbers  (the 
area  of  most  interest) . 


The  present  investigation  was  undertaken  to  study  the 
fluctuating  velocity  field  which  arises  in  the  developing 
region  of  a  rigid  walled  conduit  due  to  an  imposed  sinusoidal 
pressure  gradient.   Additionally/  results  were  sought  con- 
cerning any  phase  differences  between  the  velocity  and 
pressure  signals  with  a  view  toward  later  stability 
determinati  ons . 


II.   THEORETICAL  ANALYSIS 

A.   SOLUTION  #  1  OF  PULSATING  FLOW   EQUATIONS 

For  oscillatory   flow  in  a  circular  pipe,  the  Navier- 
Stokes  equations  must  be  solved  for  u  and  v,  the  fluctuating 
velocity  components,  as  functions  of  three  variables  (r,x,t). 
The  velocities,  u  and  v,  are  independent  of  cj> ,  the  azimuthal 
coordinate.   The  applicable  equations  from  Schlichting  [6] 
are : 

du  3u      du  1  dP    Cl  \d2u         1  du         92r' 

dt       dX       dr       p  dX 


2  2 

dr     r  8r    9x 


9v      3v      9v      1  9P 

dt      9x      dr      p  8x      J  a  2 


The  continuity  equation  is: 

It  is  assumed  that  u=u+u,  p=p+P/  and  v  =  v  where 
u  and  p  are  the  steady  components  and  u,  v,  and  p  are  the 
pulsating  components. 

After  linearization,  two  additional  assumptions  are  made, 
namely,  (1)  that  p  is  applied  across  the  full  cross-section 
of  the  tube,  i.e.,  p  is  independent  of  r ,   and  (2)  that  the 
gradient  of  v  in  the  x  direction  is  mild  compared  with  the 
gradient  in  the  r  direction. 
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The  resulting  equations  are: 


3v    -  3v 
3t      3x 


\^i  + 


V 


3r" 


r  3r 


(4) 


3u    -  3u    ~  3u 

3t      3x      3r 


p  3x 


(5) 


e*s*f-« 


(6) 


Substitution  of  (6)  into  (5)  yields 


n2  ~    1  3P    3u 

VV    U-—  7T—  -Trr 

p  3x    3t 


-  3v   uv  ,  ~  3u 

-  U  rr— +  V  7T— 

3r    r     .  3r 


(7) 


The  left  side  of  equation  (7)  is  recognized  as  that  treated 
by  Uchida  [2]  and  the  right  side  of  (7)  is  known  if  v  and  u 
are  known. 

By  assuming  v  =  R(r)X(x)e      ^    ,    equation  (4)  may  be 
solved  by  separation  of  variables.   This  technique  yields 


X_ 
X 


V 

u 


R' 


R 


+ 


RT 

rR 


=  -  w 


(8) 


for  which  the  X  solution  is 


X  =  C  e 


2 

-co  x 


(9) 


and  the  R  solution  after  assuming  a  series  representation 
becomes 


jr2    _    (sa-B2!         4  3      [ 2a  2 . 

aor    I1        —        ll92~i      r      +   192    |T"        Mr6 


|.«£-   Sli  +   tl 


192         10         48  80 


(10) 
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Therefore  the  complete  solution  for  v  is 


-oo  x   1  Nt+cM 
v  =  C,  a  r  e      e 
1  o 


3r2     8a-32  I  4 
192  Jr 


a2   a82     B3  I 
192    r  "  8     r   +  T92    10-48-   +  80  '  r 


where 


-co  U 
a  =  /  (12) 

vr 

2 

a)  U  -IN 

3  =  % (13) 

2 

and  -co   is  the  constant  introduced  in  the  separation  of 

variable  technique.   Now  that  a  solution  for  v  is  known, 
equation  (7)  may  be  solved.   By  substituting  equation  (11) 
into  (7)  and  further  assuming 

u  =  b   +  b.r  +  b0  r2  (14) 

o     1      2 


1    9p                ^         int      ~         ~,    \  int                                                  /-ir\ 

-—  •r£1=-vPne         ,u   =    u(r)  e                                                       (15) 

p  dX 

equation  (7)  becomes  a  zero  order  non-homogeneous  Bessel 
equation  of  the  form 

U"  +  i  S'.iS  u  =  PM  +  ^-^  v  |H  -  u  -1  -  Hi  >      (16) 

r     v       N    v  ^r      r    r 

The  relationship  in  equation  (14)  allows  determinations  in 

the  developing  region  to  be  made.  The  relationships  in 
equation  (15)  concerns  an  oscillating  axial  pressure  gradient 

and  a  resultant  oscillatory  axial  velocity  fluctuation. 
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The  complimentary  solution  of  equation  (16)  is 

Uc  =  C2J0(z)  +  C3Y0(z) 

where  z  =  /^T  Ar  and  X  =  J  —    . 

The  particular  solution  is  obtained  by  the  method  of 
variation  of  parameters  (Wylie[7]). 


(17) 


Up  =  Ulyl  +  U2y2 


(18) 


whe 


re  y1  =  JQ(z)  and  y2   =  Yo  ^  '  and 


-y. 


u. 


yxy2  -y2yx 


R(r,x,t) 


(19) 


u. 


y-^2  -y2yx 


R(r,x,t) 


(20) 


with 


R(r,x,t)  =  PN  + 


-int 
e     /~  9u   -  9v   uv 

v  97  "  u  8F  -  F 
v 


(21) 


The  denominator  of  u  '  and  xx'    is  the  Wronskian  of  (y, /y2) 
(Abramowitz  and  Stegun  [8]) 


W(yi*y2>  =  k 


(22) 


Therefore 


u. 


-7TZ 


2    0 


Yn(z)  R(r,x,t) 


(23) 


u. 


7TZ 


J0(z)  R(r,x,t) 


(24) 
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Entering  the  formulations  for  u  and  v  into  the  expression  for 
R(r,x,t)  and  then  substituting  the  resultant  expression 
equations  (23)  and  (24)  yields 

ttP, 


u,  =  e 


it 


N 


,1   Y0(z) 


(25) 


+  X1      z  YQ(z)  +  A2 


f  2 

z   YQ(2) 


"  X3   Z   VZ)  "  X4 


(    4 
z   YQ(z) 


X5  I2'  Y0(Z)  -A6   z6  Y0(Z) 


z7  Y0(z)  -  X7       z7  YQ(z) 


A9   z   YQ(z)  + 


...] 


where  the  various  constants  are  listed  in  Table  I 
Similarly 

ttP, 


id) 
u2  =  e 


N 


,z  J0(z)  - 


(26) 


z  JQ(z)  -  A2 


z  JQ  (z) 


A3  I z   JQ (z)  +  A4   z   JQ (z) 


f  5 

+  A5   z   JQ  (z)  +  A( 


z   Jq(z) 


+  A, 


z   Jq(z)  +  Ag   z   J0(z) 


+  A9   z   J0(z)  + 


...] 
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After  integration  to  obtain  the  expressions  for  u,  and  u„  and 
substitution  back  into  equation  (18) ,  the  complete  expression 
for  u  is  obtained. 


u  =  C2Vz)  +PNei*-X.ei*  | 


(27) 


A2e 


ic}> 


+  YQ(z)    JQ(z.) 


IT'    Jo(2)   Yo( 


+  X3e- 


2  /  2  ,, 

-  (z  -4) 


4 


i<J> 


2    3 

-  (z  -9z) 

IT 


+  9 


+  Xce 


b 


+  X_e 


■f  X8e 


J0(z)  - 


Y0(z) 


i* 


H 


i(j> 


i* 


2    4    2 

-  (z  -6z  -16) 

7T 


-  (z5-25z3+225z)  +  225 


-  (z6-36z4+208z2-96) 

7T 


J 


2    7     5       3 
-  (z  -49z  +1225zJ) 

7T 


ll,025z)  +  11,025 


r 


YQ(Z) 


j0(z: 


J0(z) 


Y0(z) 


+  A9e 


ic}) 


-  (z8-6  4z6+2  30  4z4-3  8,86  4z2+14  7,4  56) 


The  constant  C~  is  made  equal  to  zero  by  the  boundary 
condition  that  u  must  be  everywhere  finite.  The  constant 
Cp  is  solved  for  by  the  boundary  condition  that  u  =  0  at 

r  =  r0' 


15 


The  difficulty  involved  with  this  expression  is  obvious. 
Thus  a  second  attempt  at  a  theoretical  solution  was  made  to 
obtain  a  more  manageable  solution  for  the  ease  of  calculation, 

B.   SOLUTION  #  2  OF  PULSATING  FLOW  EQUATIONS 

In  this  second  phase  of  the  theoretical  analysis,  the 
assumption  that  v  is  independent  of  x  is  made.   Thus, 
equation  (4)  becomes 

lx  =  v  |af|  +  i  si  _  _^y  (28) 

8t      I  dr  r  3r    r' 

Again,  an  oscillatory  form  for  v  (now  independent  of  x) 

j  i       -^  /  \   i(Nt+<J>)    _  ,  ,  .  .  .  .     c    .-.   ■ 

is  assumed,  namely  v  =  R(r)  e        .   Substitution  of  this 

formulation  into  equation  (28)  yields  a  first  order  Bessel 

equation 

R"  +  i  R'  -  j ij  +  Ml 

r  \  r  v 

with   the    solution 

R(r)    =   J±    (i3/2    X    r)  (30) 

Incorporating  this  new  expression  for  v  alters  the  results 
for  u,  and  u9  of  the  particular  solution  of  u.   The  new 


expression  for  u„  becomes 


16 


u. 


2 z    Ji^z) H z      J0(z)    Ji(z) 


(31) 


id) 

+   e    Y 


i<|> 

+   e    Y 


X  A       z4         z2 

A24         A23Z  ^_ 

2  3  2  22 


A22Z' 


,       _2  4 

A21Z  .  2         A23Z  .  2 

2  A22Z  3  A23Z 


J02(z) 


A   z' 


J12(z) 


As  J  (x)  and  Y  (x)  are  both  cylinder  functions,  the 

expression  for  u1  can  similarly  be  expressed  by  substituting 

Y.  for  J..   The  constants  A„  .  are  tabulated  in  Table  II.   The 
11  2j 

final  expression  for  u  thus  becomes 


u  =  C2JQ(z)  +  u±   JQ(z)  +  u2  YQ(z) 


(32) 


The  immensity  of  u  has  been  somewhat  reduced.   However, 
it  was  felt  that  additional  assumptions  would  yield  a  very 
workable  solution. 

C.   SOLUTION  f  3  OF  PULSATING  FLOW  EQUATIONS 

After  obtaining  the  two  previous  solutions  for  u,  a  more 
simplified  expression  was  still  desired.   An  order  of  magni- 
tude analysis  for  u  near  the  wall  and  for  u  near  the  center- 
line  was  made,  these  being  the  two  main  areas  of  interest 
during  the  investigation.   These  order  of  magnitude  determin- 
ations did  lead  to  the  finding  of  a  simplified  solution 
applicable  across  the  entire  cross  section  of  the  pipe. 

The  expression  for  v  is 


_    i(Nt+tf))  _  ,  . 
v  =  C,  e        Ji   ^z) 


(33) 
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Based  on  the  boundary  condition  that  v  =  0  at  r  =  rn 

•    R  Cl  ei(Nt+*>  Jl(i3/2  Ar0)  =  0  (34) 

which  is  satisfied  either  by  C,  =0  or  by 

Rei(Nt+*»  Jl(i3/2  Ar0)  -  0  (35) 

Neglecting  the  trivial  solution,  then 

R   ei(Nt+«  Jl(i3/2  Ar0)  =  0  (36) 

which  leads  to 


cos(Nt+<f>)  berx  ArQ  -  sin(Nt+4>)  bei  XrQ      =      0         (37) 
or  simplified 


ber  (x  ) 

tan(Nt+*>  ■  Eii^r°T  (38) 


where  x0  =  Xrfl .   This  form  suggests  a  phase  angle  dependence 
upon  both  time  and  frequency.   For  t=0,  the  dependence  of 
phase  angle  upon  frequency  of  oscillation  is-  shown  in  Figure 
1.   The  constant  C,  in  this  eigenvalue  formulation  is 
indeterminate  because  of  the  linearity  of  the  problem. 
However,  as  in  all  such  cases,  it  is  the  form  rather  than  the 
amplitude  of  the  perturbation  velocity  which  is  of  interest. 
With  <j)  thus  determined  from  these  eigenvalue  considerations, 
v  is  of  the  form  v  =  e        Ji  ^z^  where  the  arbitrary 
constant  has  been  set  equal  to  unity  for  calculation  purposes 
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Determinations  upon  u  arise  from  equation  (7) 

%1  „2  -.    1  3P    3u      -  3v    uv    ~  3u  ,-ft% 

V  V    U  -  —  tt—  -  rr  =  ~  U  7-  -  +Vtt—  (39) 

p  3x    3t        3r    r       3r  v   ' 

Near  the  wall,  the  velocities  u  and  v  appearing  on  the  right 
side  of  equation  (39)  approach  zero.   The  velocity  gradients, 

■jr—  and  7j—  ,  though  large,  are  finite.   Thus,  near  the  wall, 
equation  (7)  becomes  that  treated  by  Uchida,  namely, 


v  v2  a  -  i|£  -  |s    =    o  (40) 

p  3x    3t 


The  resultant  solution  of  (40)  with  x  =  Ar  is 

mt  yjv±3/2x>     1S 

UWALL  "  e     —   7—372—  " 
lN   [ JQ (1    xQ) 

Near  the  centerline,  u  is  maximum,  v  is  zero  and  — 

^      i<t> 

le  3~ 

approaches  a  value  of   — ~ —  .   Likewise,  here  ^—  equals  zero 

3v  Clel* 

and  -r—  approaches  a  value  of  — * .   Thus,  near  the 

dr  z 

centerline,  equation  (7)  reduces  to 

_2  ~    1  3p    3u      -     _   icf)  ,„os 

p  dx    3t       MAX   1 

The  solution  of  this  non-homogeneous  Bessel  equation  is 

_    int   f  _  ,-3/2  .  _  .  icj)- 

Ve     V1  x)  ,    ,  ciie  "max  ; 
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In   this  expression  for  u,  the  last  term  is  negligible 
compared  to  the  first  two.   Thus,  a  complete  solution  for  u 
applicable  to  the  entire  pipe  cross-section  is 

u  =  -2 {    ° 

iN  [J0(ij/Zx0) 

This  last  form  is  equivalent  to  Uchida's  final  result  for 
the  axial  fluctuation  velocity  component.   Additionally,  the 
companion  normalized  (C,  =  1)  radial  fluctuation  velocity 
component  is 

i(nt+0)  J_  (i3/2x)  (45) 

v  =  e         1 
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III.   EXPERIMENTAL  ANALYSIS 

A.   EXPERIMENTAL  APPARATUS 

The  experimental  trials  were  conducted  in  a  horizontal 
Lucite  tube  1.5  inches  I.D.  and  45  feet  (300  diameters)  long. 
A  centrifugal  fan  was  used  to  pull  air  from  the  stagnation 
conditions  in  the  plenum  chamber  through  the  pipe  and 
pulsator  assembly.   The  inlet  chamber  consisted  of  a  rectang- 
ular box  three  feet  square  in  cross-section  and  four  feet  in 
length;  the  inlet  chamber  was  fitted  internally  with  five 
6-24  mesh  wire  screens  to  eliminate  inlet  turbulence  effects. 
Joining  the  plenum  chamber  and  the  tube  was  a  two  foot  long 
settling  chamber.   The  well  rounded  inlet  horn  was  a  quarter- 
section  of  a  9  inch  by  18  inch  ellipse  molded  of  fiber  glass 
and  polished  smooth.   The  tube,  rectangular  inlet  chamber, 
and  settling  chamber  were  leveled  and  aligned  to  eliminate 
any  curvature  effect. 

The  pulsator  assembly  consisted  of  a  sinusoidal  oscillating 
piston  arrangement  driven  by  a  MB  1250  Exciter.   The  pulsator 
was  located  at  the  end  of  the  Lucite  pipe  immediately  prior 
to  the  centrifugal  suction  fan.   The  piston  arrangement 
consisted  of  a  12  inch  diameter  disc  reciprocating  over  an 
11  inch  diameter  opening  to  the  atmosphere  within  the  13  inch 
cubical  outlet  chamber.   The  pulsating  pressure  signal  was 
created  by  vibrating  the  piston  in  a  direction  perpendicular 
to  the  direction  of  the  main  flow  velocity. 
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Velocity  measurements  were  performed  with  a  1248-L-T 
1.5-6  Thermo-Systems  Inc.  type  cross-wire  multi-sensor  probe 
and  associated  electronics.   The  hot  wire  probe  could  be 
inserted  at  measuring  ports  located  at  16,  20,  32,  40,  48, 
56,  and  272  diameters  downstream  of  the  entrance. 

In  all  cases,  the  hot  wire  probes  were  mounted  in  a 
transverse  micrometer  mechanism,  which  was  capable  of  placing 
the  sensor  to  within  .001  inch  of  any  desired  radial  position 
except  at  the  wall.   Probe  configuration  limited  placement 
at  the  wall  to  within  .05  0  inch  frcm  the  rigid  wall. 

B.   EXPERIMENTAL  PROCEDURE 

The  developing  mean  velocity  profile  was  monitored  at 
existing  ports  in  the  Lucite  pipe  located  at  20,  40,  56,  and 
272  diameters  downstream  of  the  entrance.   The  quantitative 
data  was  recorded  at  these  four  stations  for  Reynolds  Number 
values  (R  =  U    r/p)    of  2848  and  5398  respectively. 

The  fluctuating  pressure  differentials  for  high  and 
medium  mass  flows  were  obtained  using  BRUEL  &  KJAER  \      inch 
condenser  microphones  (Model  4136) .   The  mean  pressure 
differentials  were  obtained  using  micro  manometers. 

To  measure  the  pulsating  velocity  components,  the  cross 
wire  probe  was  installed  at  station  1,  twenty  diameters 
downstream  of  the  entrance.   Radial  distribution  of  velocity 
amplitudes  and  phase  shifts  were  taken  at  sinusoidal  pressure 
oscillations  of  15  cps  and  20  cps  and  at  high  and  medium 
mass  flow  rates.   All  time-dependent  data  were  recorded  on  a 
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two-channel  Sanborn  Model  296  recorder.   The  physical  setup 

of  the  electronic  equipment  utilized  to  measure  the  oscillating 

velocity  components  and  pressures  is  shown  in  Figure  2. 

On  the  two-channel  recorder  v  and  u  data  was  recorded 
alternatively;  the  input  pressure  signal  was  continuously 
monitored. 

C.   EXPERIMENTAL  RESULTS 

Initially,  mean  flow  velocity  data  was  evaluated  at  station 
1  where  the  oscillating  velocity  determinations  were  to  be 
obtained.   A  smooth  curve  was  fitted  to  the  mean  flow  data  by 
the  least  squares  method.   Figures  3  and  4  illustrate  the 
computed  plots  for  flow  Reynolds  Numbers  of  2848  and  5348 
respectively.   Table  III  gives  the   numerical  results  for  the 
developing  velocity  profile  at  the  succeeding  stations. 

Secondly,  verification  of  the  phase  shift,  <f>  ,  inherent  in 
the  radial  velocity  expression  was  sought.   After  choosing 
an  arbitrary  t  =  0  setting  upon  the  records  of  the  pressure 
and  velocity  signals,  a  subsequent  marking  of  parallel  times 
allowed  for  a  direct  verification  of  analytical  values  as 
shown  by  the  representative  tracings  in  Figures  13  and  14 
for  a  pulsation  frequency  of  20  cps .   While  equipment  limita- 
tions forbid  a  conclusive  verification  at  this  frequency; 
nevertheless,  the  appearance  of  the  90  degree  phase  shift 
between  the  axial  velocity  and  pressure  signals  on  the  center- 
line  which  is  evidenced  in  Figure  15  and  which  is  also 
predicted  by  theory  makes  the  radial  results  less  tentative. 
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Finally,  a  comparison  of  theoretical  and  experimental 
values  for  both  radial  and  axial  fluctuation  velocity 
components  was  determined  as  shown  in  Figures  5-12.   In  the 
case  of  the  axial  velocity  component,  the  comparison  between 
experimentally  measured  values  and  theoretically  calculated 
values  is  direct  whereas  for  the  radial  velocity  component 
the  theoretical  curves  incorporated  an  arbitrary  amplitude 
factor . 
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IV.   CONCLUSION 

The  final  theoretical  conclusions  follow  from  the  anamolous 
physical  situation  represented  by  equation  (28)  which  demon- 
strates the  viscous  decay  of  the  radial  velocity  component. 

Removing  the  x  dependence  of  v,  removed  the  production  term 

_  ^^ 
(u  ■s— )   which  was  necessary  to  generate  the  v  term  initially. 

Thus,  subsequent  theoretical  determinations  are  made  upon  a 

velocity  component  which  is  not  strictly  allowed  by  the  linear 

analysis  pursued.   This  anamoly  is  represented  by  the 

appearance  of  the  arbitrary  constant  in  the  final  expression 

for  the  radial  velocity  component. 

A  second  manifestation  of  this  anamoly  is  evident  in  the 

independence  of  time  that  exists  in  the  final  expression  for 

v  since 

(nt+cj))  =  tan   ber,  xQ/bei,  xQ 

Thus  the  initial  form  of  v  which  follows  from  this  result  must 
persist.   Both  the  form  and  persistence  (since  RMS  evaluations 
are  involved)  seem  to  be  suggested  by  the  comparison  of  the 
theoretical  and  experimental  values  for  the  radial  velocity 
component  portrayed  in  Figures  9-12.   Again  it  is  to  be  noted 
that  the  theoretical  curves  here  incorporate  an  arbitrary 
amplitude  factor.   However,  as  in  all  linear  eigenvalue 
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problems,  it  is  the  form  rather  than  the  absolute  value  of 
the  eigenfunction  which  is  important  and  which  is  verified 
experimentally  here. 

The  comparison  of  values  in  the  case  of  the  axial  velocity 
component  as  expressed  in  Figures  5-8  is  more  conclusive. 
This  data,  though  limited  in  nature  verifies  previous  (Richard- 
son [9]),  theoretical  and  experimental  determinations  and 
thus  helps  substantiate  the  parallel  radial  data  which  was 
simultaneously  recorded. 

These  mean  determinations  portray  another  characteristic 
of  pulsatory  flows,  namely,  that  the  maximum  mean  flow  does 
not  occur  at  the  center  of  the  pipe  but  in  a  region  near  the 
wall.   The  location  of  this  region  approaches  the  wall  as 
pulsation  frequency  increases.   As  evidenced  in  both  the 
radial  and  axial  data,  the  high  value  of  the  argument  in  the 
Bessel  functions  involved  in  the  velocity  expressions  makes 
the  wall  area  values  predominant.   Thus  the  central  portion 
of  the  flow  is  recognized  as  being  practically  unaffected  by 
pulsation.   Moreover,  the  value  of  axial  fluctuation  velocity 
is  a  non-zero  constant  on  the  pipe  axis  -  another  result 
verified  by  the  data  in  Figure  5-8.   In  contrast,  the  zero 
value  of  radial  velocity  fluctuation  on  the  pipe  axis  is 
evidenced  by  Figure  16  which  portrays  the  radial  velocity 
component  recorded  with  the  same  gain  readings  as  pertain  to 
all  other  tracings  shown.   The  appearance  of  this  figure  is 
that  of  a  filtered  noise  signal  and  demonstrates  the  favorable 
signal  to  noise  ratio  evidenced  by  the  other  recorded  data. 
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In  addition  to  the  results   concerning  the  existence  of 
a  radial  fluctuation  velocity  component,  the  second  new 
result  warranted  by  the  accompanying  analysis  was  the  appear- 
ance of  a  frequency  dependent  phase  shift  in  the  formulation 
of  this  velocity  component.   Results  concerning  this  phase 
shift  are  shown  in  Figures  13  and  14  where  determinations 
for  t  =  0  values  approach  those  values  of  phase  shift 
suggested  by  the  analysis  for  a  pulsation  frequency  for  20  cps. 
Equipment  limitations  did  not  allow  verification  of  non-zero 
time  values.   However,  as  stated  previously,  the  appearance 
of  a  90  degree  phase  shift  in  the  axial  velocity  and  pressure 
data  recorded  in  Figure  15  tends  to  substantiate  the  companion 
radial  velocity  data. 

In  conclusion,  the  limited  experimental  program  reported 
here  only  demonstrates  trends  toward  the  parallel  theoretical 
determinations.   A  more  comprehensive  experimental  program 
would  require  both  new  electronic  filters  with  band  pass 
capabilities  in  the  .01  cps  range  and  a  new  pressure  monitoring 
device  to  replace  the  microphone  system  which  has  severely 
attenuated  response  below  20  cps.   In  view  of  these  require- 
ments, a  re-evaluation  of  analytical  results  is  recommended 
before  a  more  extensive  experimental  phase  is  commenced. 
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TABLE  I 

Constants  for  series  solution  of  v;  solution  #  1  of  the 
Pulsating  Flow  equations. 
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TABLE  II 

Constants  for  v  independent  of  x  solution;  solution  #  2 
of  the  Pulsating  Flow  equations . 
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TABLE  III 


Coefficients  of  the  polynomial  approximations  to  the  velocity 
profiles  of  the  form: 
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\  2 


+  b3rT 


x/d   mass  flow  rate 


20 


40 


56 


272 


high 

medium 

high 

medium 

high 

medium 

high 

medium 


_P_ 
0.9915 

0.9760 

1.0137 

0.9804 

1.0010 

0.9996 

0.9653 

0.9904 


0964 
3159 
4726 
1590 
1717 
1579 
6795 
1320 


1.3958 
-.1252 
2.3368 
0.5294 
1.3127 
0.9104 
1.6906 
-1.5630 


-2.4025 
-1.2745 
-2.9733 
-1.7836 
-2.1629 
-1.8154 
-0.0106 
0.4163 


31 


20 


10 


-10 


-20 


-30 


-40 


i 

1 

! 

1 

"  "  i 

<j:u 

..j 

.:.. 

.10         15 
Frequency  (cps) 


20 
— rn3T 


Figure  1.   Phase  Shift  (degrees)  vs  Frequency  (cps)  of 
v  compared  with  pressure  signal  at  t  =  0. 
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Figure  2.   Physical  arrangement  of  electronic 
equipment  for  measuring  u,  v,  and  p 
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Figure  3.   Developing  Velocity  Profile;  x/d=20, 
R=2848. 
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Figure  4.   Developing  Velocity  Profile;  x/d=20, 
R=5348. 
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Figure  5.   Axial  Velocity  Component;  R=5348, 
Freq  =  20  cps . 
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Figure    6.      Axial   Velocity   Component;    R=5348, 
Freq   =    15    cps 
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Figure    7.      Axial   Velocity   Component;    R=2848, 
Freq   =    15    cps . 
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Figure  8.   Axial  Velocity  Component;  R=2848, 
Freq  =  15  cps . 
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Figure  9.   Radial  Velocity  Component;  R=534  8, 
Freq  =  20  cps. 
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Figure  10.   Radial  Velocity  Component;  R=534  8, 
Freq  =  15  cps . 
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Figure  11.   Radial  Velocity  Component;  R=2  84  8, 
Freq  =  20  cps . 
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Figure  12.   Radial  Velocity  Component;  R=2  84  8, 
Freq  =  15  cps . 
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Figure  13.   Recorder  Tracing  of  p  (left)  and 
v  (right)  Radius  =  0.7,  R=2848,  Freq  = 
20  cps,  4>  ~  45°. 
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Figure  14.   Recorder  Tracing  of  p  (left)  and 
v  (right)  Radius  =  0.7,  R=534  8,  Freq  = 
20  cps,  <f)  ~  45°. 
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Figure  15.   Recorder  Tracing  of  p  (left)  and 
U  (right)  Radius  =0.0,  R=5348,  Freq  = 
20  cps,  <f>  ~  90°  . 
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Figure  16.  Recorder  Tracing  p  (left)  and 
v  (right)  Radius  =0.0,  R=5348,  Freq  = 
15  cps ,  v  *  0  . 
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ABSTRACT 

A  study  of  the  fluctuating  velocity  field  in  the  developing 
region  of  a  pulsating  pipe  flow  is  presented. 

Theoretically,  three  separate  determinations  of  radial  and 
axial  perturbation  velocity  profiles  were  calculated  involving 
varying  degrees  of  analytical  simplication. 

Experimentally,  pulsating  disturbances  were  created  by  a  sinu- 
soidal pressure  wave  exciter  installed  at  the  outlet  of  a  1.5  inch 
circular  wind  tunnel  facility.  Cross-wire  anemometer  measurements 
of  the  pulsating  velocity  components  compared  favorably  with  the 
parallel  theoretical  results.  Likewise,  experimental  determinations 
of  the  phase  angles  between  velocity  and  pressure  signals  compared 
favorably  with  analytical  values. 
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